Comparison of machine learning algorithms used to classify the asteroids observed by all-sky surveys

Context. Multifilter photometry from large sky surveys is commonly used to assign asteroid taxonomic types and study various problems in planetary science. To maximize the science output of those surveys, it is important to use methods that best link the spectro-photometric measurements to asteroid taxonomy. Aims. We aim to determine which machine learning methods are the most suitable for the taxonomic classification for various sky surveys. Methods. We utilized five machine learning supervised classifiers: logistic regression, naive Bayes, support vector machines (SVMs), gradient boosting, and MultiLayer Perceptrons (MLPs). Those methods were found to reproduce the Bus-DeMeo taxonomy at various rates depending on the set of filters used by each survey. We report several evaluation metrics for a comprehensive comparison (prediction accuracy, balanced accuracy, F1 score, and the Matthews correlation coefficient) for 11 surveys and space missions. Results. Among the methods analyzed, multilayer perception and gradient boosting achieved the highest accuracy and naive Bayes achieved the lowest accuracy in taxonomic prediction across all surveys. We found that selecting the right machine learning algorithm can improve the success rate by a factor of >2. The best balanced accuracy (similar to 85% for a taxonomic type prediction) was found for the Visible and Infrared Survey telescope for Astronomy (VISTA) and the ESA Euclid mission surveys where broadband filters best map the 1 mu m and 2 mu m olivine and pyroxene absorption bands. Conclusions. To achieve the highest accuracy in the taxonomic type prediction based on multifilter photometric measurements, we recommend the use of gradient boosting and MLP optimized for each survey. This can improve the overall success rate even when compared with naive Bayes. A merger of different datasets can further boost the prediction accuracy. For the combination of the Legacy Survey of Space and Time and VISTA survey, we achieved 90% for the taxonomic type prediction.Peer reviewe

Near infra-red spectroscopy of the asteroid 21 Lutetia. II. Rotationally resolved spectroscopy of the surface

Reproduced with permission. Copyright ESO. Article published by EDP Sciences and available at http://www.aanda.orgInternational audienceAims. In the framework of the ground-based science campaign dedicated to the encounter with the Rosetta spacecraft, the mineralogy of the asteroid (21) Lutetia was investigated. Methods. Near-infrared (NIR) spectra of the asteroid in the 0.8−2.5 μm spectral range were obtained with SpeX/IRTF in remote observing mode from Meudon, France in March and April 2006. We analysed these data together with previously acquired spectra - March 2003, August 2004. I-band relative photometric data obtained on 20 January 2006 using the 105 cm telescope from Pic du Midi, France has been used to build the ephemeris for physical observations. A χ2 test using meteorite spectra from the RELAB database was performed in order to find the best fit of complete visible + infrared (VNIR) spectra of Lutetia. Results. The new spectra reveal no absorption features. We find a clear spectral variation (slope), and a good correspondence between spectral variations and rotational phase. Two of the most different spectra correspond to two opposite sides of the asteroid (sub-Earth longitude difference around 180◦). For the neutral spectra a carbonaceous chondrite spectrum yields the best fit, while for those with a slightly positive slope the enstatitic chondrite spectra are the best analog. Based on the chosen subset of the meteorite samples, our analysis suggests a primitive, chondritic nature for (21) Lutetia. Differences in spectra are interpreted in terms of the coexistence of several lithologies on the surface where the aqueous alteration played an important role

The non-convex shape of (234) Barbara, the first Barbarian

Asteroid (234) Barbara is the prototype of a category of asteroids that has been shown to be extremely rich in refractory inclusions, the oldest material ever found in the Solar System. It exhibits several peculiar features, most notably its polarimetric behavior. In recent years other objects sharing the same property (collectively known as "Barbarians") have been discovered. Interferometric observations in the mid-infrared with the ESO VLTI suggested that (234) Barbara might have a bi-lobated shape or even a large companion satellite. We use a large set of 57 optical lightcurves acquired between 1979 and 2014, together with the timings of two stellar occultations in 2009, to determine the rotation period, spin-vector coordinates, and 3-D shape of (234) Barbara, using two different shape reconstruction algorithms. By using the lightcurves combined to the results obtained from stellar occultations, we are able to show that the shape of (234) Barbara exhibits large concave areas. Possible links of the shape to the polarimetric properties and the object evolution are discussed. We also show that VLTI data can be modeled without the presence of a satellite.Comment: 10 pages, 6 figure

(16) Psyche: A mesosiderite-like asteroid?

Asteroid (16) Psyche is the target of the NASA Psyche mission. It is considered one of the few main-belt bodies that could be an exposed proto-planetary metallic core and that would thus be related to iron meteorites. Such an association is however challenged by both its near- and mid-infrared spectral properties and the reported estimates of its density. Here, we aim to refine the density of (16) Psyche to set further constraints on its bulk composition and determine its potential meteoritic analog. We observed (16) Psyche with ESO VLT/SPHERE/ZIMPOL as part of our large program (ID 199.C-0074). We used the high angular resolution of these observations to refine Psyche's three-dimensional (3D) shape model and subsequently its density when combined with the most recent mass estimates. In addition, we searched for potential companions around the asteroid. We derived a bulk density of 3.99\,$\pm$\,0.26\,g$\cdot$cm$^{-3}$ for Psyche. While such density is incompatible at the 3-sigma level with any iron meteorites ($\sim$7.8\,g$\cdot$cm$^{-3}$), it appears fully consistent with that of stony-iron meteorites such as mesosiderites (density $\sim$4.25\,$\cdot$cm$^{-3}$). In addition, we found no satellite in our images and set an upper limit on the diameter of any non-detected satellite of 1460\,$\pm$\,200}\,m at 150\,km from Psyche (0.2\%\,$\times$\,R$_{Hill}$, the Hill radius) and 800\,$\pm$\,200\,m at 2,000\,km (3\%\,$\times$\,$R_{Hill}$). Considering that the visible and near-infrared spectral properties of mesosiderites are similar to those of Psyche, there is merit to a long-published initial hypothesis that Psyche could be a plausible candidate parent body for mesosiderites.Comment: 16 page

Physical properties of ESA Rosetta target asteroid (21) Lutetia: Shape and flyby geometry

Aims. We determine the physical properties (spin state and shape) of asteroid (21) Lutetia, target of the ESA Rosetta mission, to help in preparing for observations during the flyby on 2010 July 10 by predicting the orientation of Lutetia as seen from Rosetta. Methods. We use our novel KOALA inversion algorithm to determine the physical properties of asteroids from a combination of optical lightcurves, disk-resolved images, and stellar occultations, although the latter are not available for (21) Lutetia. Results. We find the spin axis of (21) Lutetia to lie within 5 degrees of ({\lambda} = 52 deg., {\beta} = -6 deg.) in Ecliptic J2000 reference frame (equatorial {\alpha} = 52 deg., {\delta} = +12 deg.), and determine an improved sidereal period of 8.168 270 \pm 0.000 001 h. This pole solution implies the southern hemisphere of Lutetia will be in "seasonal" shadow at the time of the flyby. The apparent cross-section of Lutetia is triangular as seen "pole-on" and more rectangular as seen "equator-on". The best-fit model suggests the presence of several concavities. The largest of these is close to the north pole and may be associated with large impacts.Comment: 17 pages, 5 figures, 3 tables, submitted to Astronomy and Astrophysic

Physical and dynamical properties of the unusual V-type asteroid (2579) Spartacus

Context. Asteroid (2579) Spartacus is a small V-type object located in the inner main belt. This object shows spectral characteristics unusual for typical Vestoids, which may indicate an origin deeper than average within Vesta or an origin from an altogether different parent body. Aims. Our main goal is to study the origin of Spartacus. We derive the spin of Spartacus and a convex shape model of Spartacus in order to increase the knowledge of the body's physical properties. The rotational parameters are then used to investigate dynamical evolution of the object as well as to distinguish regions sampled by spectral observations to determine whether its surface displays heterogeneity. Methods. We collected lightcurves available from the literature (oppositions of 2009, 2012) and obtained additional photometric observations at various telescopes in 2016, 2017, and 2018. We used the lightcurve inversion method to derive a spin and convex shape model. We have collected spectral observations over two rotational periods of Spartacus and determined its spectral parameters using the modified Gaussian model (MGM). We then dynamically integrated the orbital elements of Spartacus, taking into account existing information, including its thermal properties, size and the derived spin axis orientation. Results. We find two models for (2579) Spartacus: (a) lambda = 312 degrees +/- 5 degrees, beta = -57 degrees +/- 5 degrees and (b) lambda = 113 degrees +/- 5 degrees, beta = -60 degrees +/- 5 degrees both retrograde. We find that the drift direction for Spartacus is consistent with separation from Vesta, and after a backward integration of 1 Gyr the asteroid reaches the boundary of the family. We did not observe spectral variations with rotation, thus the body most likely has a homogeneous surface. Additionally, new spectral analysis indicates that the 1.0 and 2.0 mu m band centers are within ranges that are typical for Vestoids while the area ratio of these bands is about half that of typical Vestoids. Conclusions. The asteroid (2579) Spartacus is in retrograde rotation and has a drift direction consistent with an origin from Vesta. The revised spectral band centers are within ranges typical for Vestoids, while band area ratio (BAR) is unusually low compared to that of other V-types. The dynamical model shows that the asteroid could have migrated to its current location from the edges of the Vesta family within 1 Gyr, but an origin from an earlier impact on Vesta could also be plausible.Peer reviewe

Radar Observations and the Shape of Near-Earth Asteroid 2008 EV5

We observed the near-Earth asteroid 2008 EV5 with the Arecibo and Goldstone planetary radars and the Very Long Baseline Array during December 2008. EV5 rotates retrograde and its overall shape is a 400 /pm 50 m oblate spheroid. The most prominent surface feature is a ridge parallel to the asteroid's equator that is broken by a concavity 150 m in diameter. Otherwise the asteroid's surface is notably smooth on decameter scales. EV5's radar and optical albedos are consistent with either rocky or stony-iron composition. The equatorial ridge is similar to structure seen on the rubble-pile near-Earth asteroid (66391) 1999 KW4 and is consistent with YORP spin-up reconfiguring the asteroid in the past. We interpret the concavity as an impact crater. Shaking during the impact and later regolith redistribution may have erased smaller features, explaining the general lack of decameter-scale surface structure.Comment: This paper has been accepted for publication in Icarus: http://www.sciencedirect.com/science/article/B6WGF-5207B2F-4/2/d87cd2ae4da00c2b277e2dc79a532c4

(704) Interamnia: a transitional object between a dwarf planet and a typical irregular-shaped minor body

Context. With an estimated diameter in the 320–350 km range, (704) Interamnia is the fifth largest main belt asteroid and one of the few bodies that fills the gap in size between the four largest bodies with D > 400 km (Ceres, Vesta, Pallas and Hygiea) and the numerous smaller bodies with diameter ≤200 km. However, despite its large size, little is known about the shape and spin state of Interamnia and, therefore, about its bulk composition and past collisional evolution. Aims. We aimed to test at what size and mass the shape of a small body departs from a nearly ellipsoidal equilibrium shape (as observed in the case of the four largest asteroids) to an irregular shape as routinely observed in the case of smaller (D ≤ 200 km) bodies. Methods. We observed Interamnia as part of our ESO VLT/SPHERE large program (ID: 199.C-0074) at thirteen different epochs. In addition, several new optical lightcurves were recorded. These data, along with stellar occultation data from the literature, were fed to the All-Data Asteroid Modeling algorithm to reconstruct the 3D-shape model of Interamnia and to determine its spin state. Results. Interamnia’s volume-equivalent diameter of 332 ± 6 km implies a bulk density of ρ = 1.98 ± 0.68 g cm−3, which suggests that Interamnia – like Ceres and Hygiea – contains a high fraction of water ice, consistent with the paucity of apparent craters. Our observations reveal a shape that can be well approximated by an ellipsoid, and that is compatible with a fluid hydrostatic equilibrium at the 2σ level. Conclusions. The rather regular shape of Interamnia implies that the size and mass limit, under which the shapes of minor bodies with a high amount of water ice in the subsurface become irregular, has to be searched among smaller (D ≤ 300 km) less massive (m ≤ 3 × 1019 kg) bodies

(216) Kleopatra, a low density critically rotating M-type asteroid

Context. The recent estimates of the 3D shape of the M/Xe-type triple asteroid system (216) Kleopatra indicated a density of ~5 g cm−3, which is by far the highest for a small Solar System body. Such a high density implies a high metal content as well as a low porosity which is not easy to reconcile with its peculiar “dumbbell” shape. Aims. Given the unprecedented angular resolution of the VLT/SPHERE/ZIMPOL camera, here, we aim to constrain the mass (via the characterization of the orbits of the moons) and the shape of (216) Kleopatra with high accuracy, hence its density. Methods. We combined our new VLT/SPHERE observations of (216) Kleopatra recorded during two apparitions in 2017 and 2018 with archival data from the W. M. Keck Observatory, as well as lightcurve, occultation, and delay-Doppler images, to derive a model of its 3D shape using two different algorithms (ADAM, MPCD). Furthermore, an N-body dynamical model allowed us to retrieve the orbital elements of the two moons as explained in the accompanying paper. Results. The shape of (216) Kleopatra is very close to an equilibrium dumbbell figure with two lobes and a thick neck. Its volume equivalent diameter (118.75 ± 1.40) km and mass (2.97 ± 0.32) × 1018 kg (i.e., 56% lower than previously reported) imply a bulk density of (3.38 ± 0.50) g cm−3. Such a low density for a supposedly metal-rich body indicates a substantial porosity within the primary. This porous structure along with its near equilibrium shape is compatible with a formation scenario including a giant impact followed by reaccumulation. (216) Kleopatra’s current rotation period and dumbbell shape imply that it is in a critically rotating state. The low effective gravity along the equator of the body, together with the equatorial orbits of the moons and possibly rubble-pile structure, opens the possibility that the moons formed via mass shedding. Conclusions. (216) Kleopatra is a puzzling multiple system due to the unique characteristics of the primary. This system certainly deserves particular attention in the future, with the Extremely Large Telescopes and possibly a dedicated space mission, to decipher its entire formation history

Evidence for differentiation of the most primitive small bodies

Context. Dynamical models of Solar System evolution have suggested that the so-called P- and D-type volatile-rich asteroids formed in the outer Solar System beyond Neptune's orbit and may be genetically related to the Jupiter Trojans, comets, and small Kuiper belt objects (KBOs). Indeed, the spectral properties of P- and D-type asteroids resemble that of anhydrous cometary dust. Aims. We aim to gain insights into the above classes of bodies by characterizing the internal structure of a large P- and D-type asteroid. Methods. We report high-angular-resolution imaging observations of the P-type asteroid (87) Sylvia with the Very Large Telescope Spectro-Polarimetric High-contrast Exoplanet REsearch (SPHERE) instrument. These images were used to reconstruct the 3D shape of Sylvia. Our images together with those obtained in the past with large ground-based telescopes were used to study the dynamics of its two satellites. We also modeled Sylvia's thermal evolution. Results. The shape of Sylvia appears flattened and elongated (a/b 1.45; a/c 1.84). We derive a volume-equivalent diameter of 271 ± 5 km and a low density of 1378 ± 45 kg m-3. The two satellites orbit Sylvia on circular, equatorial orbits. The oblateness of Sylvia should imply a detectable nodal precession which contrasts with the fully-Keplerian dynamics of its two satellites. This reveals an inhomogeneous internal structure, suggesting that Sylvia is differentiated. Conclusions. Sylvia's low density and differentiated interior can be explained by partial melting and mass redistribution through water percolation. The outer shell should be composed of material similar to interplanetary dust particles (IDPs) and the core should be similar to aqueously altered IDPs or carbonaceous chondrite meteorites such as the Tagish Lake meteorite. Numerical simulations of the thermal evolution of Sylvia show that for a body of such a size, partial melting was unavoidable due to the decay of long-lived radionuclides. In addition, we show that bodies as small as 130-150 km in diameter should have followed a similar thermal evolution, while smaller objects, such as comets and the KBO Arrokoth, must have remained pristine, which is in agreement with in situ observations of these bodies. NASA Lucy mission target (617) Patroclus (diameter ≈140 km) may, however, be differentiated